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The correlation radiometer is analyzed to determine the sensitivity that can be obtained under 
various operating conditions. 

The radiometer using a sine wave comparison signal is analyzed and compared with the usual 
radiometer that employs a random noise for the comparison signal. It is found that the radiometer 
employing the sine wave comparison signal is the more sensitive of the two circuits, particularly in 
the case that the effective temperature of the input noise signal is greater than the effective input 
temperature of the amplifiers. 

It is shown that if nonidentical amplifiers are used in the correlation circuit, the properties of the 
radiometer are determined by the portion of the amplifier response functions in the frequency interval 
that the two response functions overlap. The effect of amplifier gain fluctuations are considered, and 
although the correlation scheme reduces the effect of gain fluctuations, it is shown that they still do 
contribute to the output fluctuations of the radiometer. 

Calculations are included showing that the effect of a differential phase shift between the two chan- 
nels is a reduction in radiometer sensitivity. The same conclusion is reached concerning the effect of a 
differential time delay. 

Finally, it is shown that if the comparison signal and the input signal have the same statistical 
properties, the requirements on the multiplier are less stringent than if the two signals have different 
statistical properties. 

Key Words: Correlation radiometer, differential time delay, gain fluctuations, imperfect multi- 
plier, noise comparison, nonidentical amplifiers, sine-wave comparison signal. 



1 . Introduction 

The use of correlation techniques in radiometry has been suggested recently by several workers 
[Strum, 1958; Blum, 1959; Colvin, 1961; Allred, 1962]. One of the major reasons for this interest 
is that by using correlation techniques, it is possible to build a radiometer in which amplifier gain 
fluctuations contribute less to the fluctuations present in the output of the instrument than they do 
in the conventional radiometer; thus, the correlation radiometer might be expected to be superior 
to the conventional radiometer under conditions in which gain fluctuations are important. Strum 
[1958J has pointed out that as the noise level in a conventional radiometer system is reduced, the 
fluctuations due to changes in amplifier gain become more serious; therefore, at cryogenic tempera- 
tures it is very important to reduce this contribution to the total radiometer fluctuation. Thus, 
it was felt to be desirable to calculate the sensitivity obtainable from instruments of this type. 

The analysis was carried out on the correlation receiver shown in figure 1. Although several 
different types of correlation receivers have been suggested in the literature, the operation of all 
of them is quite similar to the one chosen here. Therefore, the results obtained here should be 
applicable with appropriate modifications to many of the correlation radiometer circuits in the 
literature. 

The first problem to be considered (see sec. 2) is the calculation of the sensitivity of a radiom- 
eter with amplifiers whose gain and passband are identical but with different effective tempera- 
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FIGURE 1. Radiometer circuit. 



tures at the input circuit. These calculations are carried out for two types of input signals: 
(1) both input signals consist of white noise and, (2) one input signal is white noise and the other 
is sinusoidal. The first case was selected because it corresponds to the problem considered by 
Colvin [1961], and the second was chosen because it was the problem considered by Allred [1962]. 
Throughout this calculation, all filters are assumed to have a square passband. This assumption 
is made because it allows the results of the calculation to be written rather explicitly and thus it 
is easy to see the physical significance of the results. When it is desired to apply these results 
to a real radiometer with filters that are not square, it is only necessary to use the formalism of 
Colvin [1959] that expresses the effective bandwidth of a filter in terms of a convolution integral. 
Thus, this generalization can be carried out rather simply. 

Section 3 is an analysis of the correlation radiometer in which it is assumed that the two ampli- 
fiers have different gains and different response curves. 

Section 4 deals with the same circuit, but here the effect of variation in phase and in time delay 
of the signals in the two amplifiers is considered. 

Finally, section 5 deals with some effects produced by an imperfect multiplier. 

2. Sensitivity of an Ideal Correlation Radiometer 

The circuit to be analyzed is shown in figure 1. Two noise sources, designated X and F, are 
connected to opposite arms of a matched hybrid junction. The signals from the other pair of arms 
are applied to amplifiers whose transfer functions are Ri((o) and Rzico) respectively. The outputs 
of the amplifiers are multiplied together and the product is filtered with a smoothing filter (low 
pass filter). The output of the smoothing filter is displayed on a d-c instrument. 

As the following analysis will show, the average value of the deflection of the output instru- 
ment is proportional to the difference in noise temperature of the two sources, X and Y. Thus, 
when the two noise sources have the same effective temperature, the average value of the output 
deflection is zero and the radiometer is said to be balanced. The instantaneous deflection is a 
stochastic function of the time; thus the random fluctuations of the output deflection produce an 
uncertainty in the experimental conditions that correspond to a balance. 

The sensitivity of a radiometer is usually defined as the change in temperature of one of the 
noise sources that will produce a deflection whose magnitude is equal to the root mean square of 
the output fluctuation. Therefore, in order to calculate the sensitivity of a radiometer it is neces- 
sary to obtain expressions for both the rms value of the output fluctuation and the change in aver- 
age deflection per degree change in effective temperature of one of the noise sources. 

In this section the assumption will be made that the two amplifiers have the same power trans- 
fer function, and for simplicity, a "square passband" will be assumed. 

If x(t) is the signal voltage due to source X that leaves each amplifier, and y(t) is the signal 
voltage due to source Fthat leaves amplifier 1, the signal voltage —y(t) due to source Fmust leave 
amplifier 2. It will be further assumed that the amplifiers introduce noise voltages Z\ and z>i 
respectively. 

Since the random functions x{t\ y(t), Z\{t\ and z 2 (t) represent noise from physically separate 
noise generators, it will be assumed that they are mutually independent. It is also assumed that 
each is a second order, stationary Gaussian process possessing a continuous spectral density 
function. These assumptions imply that the processes are ergodic so that time averages and 
ensemble averages coincide. It is also assumed that each of these functions has zero mean. 
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If the output of amplifier 1 is designated by v\ and the output of amplifier 2 by i; 2 , then 

v\=x + y-f z\ (1) 

and 

v 2 = x-y+z 2 . (la) 

The multiplier will be assumed to have a response law such that 

u = av\Vz (2) 

where u is the voltage output and a is a constant of proportionality. 

For any random quantity £, let £ be the expectation or ensemble average of f . (This is the 
equivalent of the notation E(Q used in the literature of mathematical statistics.) Using the inde- 
pendence and normality of the random processes and the fact that each process has zero mean, 
the autocorrelation function of u is 



i// w (t) = u(t )u(t + t)= a 2 vi(t)v>(t )vi(t + r)v 2 (t + r) 



Wr) - a 2 [i//,.,(r) + Wt) - 2i//,.(0)i//,/0) + x\f Zx (r) (i//,.(r) + Wt))+^ 2 (t)(Wt) + i^t)) + iM*)*rXr)]. (3) 
If fj, /= 1, 2, 3, 4 are Gaussian variables, each with zero mean and 

then 



£ 1&&& = Or 12 CT; i 4 H" CT130-24 + 0"l40*23. (4) 

Setting fi = f2 = ^) and f 3 = f* = *(* + t), there results 



Wr) - x\t)x\t + t) = 2i//|(r) + i//|(0). (5) 

When (4) is applied to (3), the result is 

a-*U?) = 2 [#(t) + *// 2 (t)] + [ifjAO) - %(0)] 2 + <K(t) [ Wr) + «r)] 

+ Wr) [Wt) + i//,(r)] + ^(t>Mt). (6) 

The Fourier transform of the autocorrelation function of a process is the spectral density 
of the process. Further, the Fourier transform of a product of functions is the convolution of the 
Fourier transforms. Thus, designating the spectral density by Q, and the convolution operation 
by *, 

a- 2 Q„(f) = [<M0) - «0)] 2 8(/) + 2 [Q x * QAf) + Qy * <?y(/)] 

+ W.t+ o J *■[(?*+ W c^+o.;* <?**(/). (?) 
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Case I. Consider both X and Y to be white noise generators at temperature T x and T y respec- 
tively. When the spectral densities are evaluated, it is necessary to consider the power division 
that takes place at the hybrid junction. These spectral densities are (see appendix) 

Q x *Q x ~±(kT x fB, (8) 

Q x *Q Zl ~\k 2 T x T Zl B, (8a) 

Qz^Qz^\k 2 T Zl T z JB. (8b) 

Thus, using Q' to designate the portion of the spectral density that contributes to fluctuations, 
this portion of (7) becomes 

Q'u = ^p [?1+ T$ + T x (T Zl + T i2 ) + T y (T Zx + T H ) + 2T Zl T Ss ] . (9) 

With the assumption that the two amplifiers have the same effective noise temperature, this 
becomes 

<?i=pp m+TI + 2T x Tz+2T u T z + 2Tn (10) 

If the smoothing filter has a power response G(/), the output power spectrum from the smoothing 

filter is 

W(f) = Cif)Q' u if). (11) 

It is well known that integrating a power spectrum over all frequencies results in the mean square 
of the output voltage or current; so if the output voltage from the smoothing filter is wit), and 
G{f) represents a square pass filter of width 6, then 

= J" W(f)df= \ a 2 k 2 G bB [71 + T\ + 27* + 2(T X T Z + T,T 2 )] (12) 



w{tf 



where Go is the power response of the filter at zero frequency. The fluctuations at balance are 
required, so the balance condition is substituted into w 2 . The balance occurs when T x ~T y , so 
for balance 

~w~ 2 = a 2 k 2 GM{T x + T z f. (13) 

Therefore, in terms of an rms voltage 

wrms = W 1/2 = ak(GobByi*(T x + T,). (14) 

Next, we calculate the deflection arising from a certain imbalance of the input signals. From (2) 

a = (IV1V2 = a(x H- y H- zi) (x — y -f 22), 

= a[x 2 — y 2 + ziz 2 + zi(x — y) + z 2 (x -f- y)]. (15) 

Since x(t), y{t), z\{t) and zz(t) are all uncorrelated, when (15) is averaged over time, the result is 

u = a(x~ 2 -f). (16) 
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Again applying the theorem that the integral of the power spectrum over all frequencies is the 
average square of a function, (16) becomes 



u = a 



JQxdf-JQydf, 



= ^f[T x -T y ]. 

In order to examine the effect of a small imbalance, set 

Ty=T X +&T. 



Then 



Since u = when T x = T y , 



u + Au=?ML[T x -(T x +M)]. 



|Au|=^A7\ 



The resulting signal at the output of the smoothing filter is 

|Aw|=CJ/ 1 |Att|=|c^*BAT. 



(17) 



(18) 



(19) 



(20) 



As is usual in radiometer calculations, the assumption is made that the minimum detectable signal 
occurs when the deflection is equal to w rms . Thus, equating the right-hand sides of (14) and (20), 
the condition corresponding to minimum detectable signal results; i.e., 



ak(G Q bBY l2 (T x + T M ) = - GfkB&F. 



Thus, the minimum detectable temperature difference is 

/\T= 2(1>I BY l 2 (T. r + T z ). 
Case II. Let Y produce a sinusoidal signal, so that 



(21) 



(22) 



Qy(f) = jW-fo) + 8(f+M. 



(23) 



The expression for ip X 2-y* is needed, so it will be calculated first. From the previous calcu- 
lation, it is known that 



$X2-y2 = *\>X2 + l/>l/2 — 2x 2 y 2 . 



Thus, it is necessary to compute i// v2 : 



if'i^rHOrV+T), 



— — COS 2 (Oot cos 2 coo(£ + t), 
4 



\\iy2-- 



3 2 ■ ! • 2 

- cos z cooT-h- snr co r 
o o 



(24) 
(25) 
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Making use of the above and (5), the result is 



— 2 v 4 /3 1 \ — v 2 

tyx2-y2 = x 2 + 2(i// x ) 2 + — - cos 2 o) T + - sin 2 o) T ) — 2x 2 — 



4 \8 



(26) 



From this, the power spectral density is 
Qx2-y2= j ip X 2-y2e- i(OT dr, 



— 2 _ ,,2 



* 2 f)o(/)+2(^*£y+^ 



\ S(/) + ^ 8(/"-2/ ) + ^ 8(/-+2/o) 



Therefore 



Qxi-tfi + Qz,(x-y) + Qz 2 (x 



+V)+QziZ2 



(27) 



(28) 



" jf S( ^ + ^ [8(/ ~ 2/o) + 8C/+ 2/o)] + 2(<?J - * ^ 



+ Qz 1 *(Q*+Qy)+ Q Z2 *(Q X +Q V )+ Q*, * <?*,} ■ 



(29) 



Since the present calculation is to compute the fluctuations occurring in this radiometer, only 
the terms that contribute to the fluctuation need be considered; therefore they will be expressed as 



& = a 2 {2(Q x *Q x ) + «?,, + QzMQ* + Qv) + <?«>*<?*}, 



(30) 



where Quif) denotes the fluctuating portion of Q u . 

All of the evaluations of the convolutions that occur in (30) can be taken from (8), except 
Qy*Qz x \ which is, 

Qv*Qz, =&,*<?«=/ Qz>(f')Qy(J-f'W' , 

where Q y is given in (23). Thus, 



Q 



v*Qz 1 =j 



Q4f) j Wf-f +/o) +«(/■-/' -/oM\ 



v' 



(31) 



=7K?^+/o)+<?.,V-/i)]. 



Since Q Zl is assumed to be constant through the passband, when/ is in the passband and /is 
small, this convolution becomes 



Qy*Qz t =J kT Zl . 



(32) 



Therefore 



« = a 2 l\ (kT x f + \ VBTJJ , , l + TJ+\ *KT tl + TJ+\ VTJjV 
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By analogy with the previous calculation, the spectral density of the output of the filter can be 
written as 

W(f) = G(j)Q u (f), 

and 



1 = j " W{f)df= <' 2 ('o(2b) {I {kT x fB + \ k*BTJT„ + TJ + 1 v 2 k(T z , + T Z2 ) + 1 VT H TJ 



(34) 



The d-c term of (29) indicates that the system is balanced when kT x B = v' 2 /2. Also, to simplify 
the expression, it will be assumed that T z , = T Z2 . Then, at balance 

J?= a 2 G n bk 2 B Q 71+ 2r. r r,+ T*j • (35) 

By analogy with the previous calculation, a root mean square fluctuation amplitude may be de- 
fined by 

w Tms =ak(GobB) i l* (| T%+2T X T,+ TJ) 1 *> (36) 

As mentioned previously, in order to arrive at the radiometer sensitivity, two parameters 
must be computed. The first is w rms , which is given in (36). The second is the change in deflec- 
tion due to a small change in temperature of the thermal source. The deflection sensitivity will 
be the same as it was in the previous calculation, so it is 

\Aw\ = G y^AT- 

Again, the minimum change in temperature that is observable will be taken to be a change 
that produces a deflection equal to the w rms . Therefore 



G i/2 ^ AT=ak(G {) Byn Q ri+2T x T z +T*\ 



1/2 



and the final result is 



AT=2 (|)" 2 (| n+2T x T z + li) U2 - (37) 

When this result is compared with (22), the sensitivity of a radiometer with both signals as- 
sumed to be white noise, it is evident that this instrument has a somewhat higher sensitivity. 
This is particularly true in the case that the input signal is the dominant noise in the system and 
thus sets the fluctuation level. In this case the improvement in sensitivity approaches V2. 

3. Radiometer With Dissimilar Amplifiers 

In section 2 the assumption was made that the two amplifiers were identical, except that they 
were allowed to have different effective temperatures. Here the more general problem, in which 
the two amplifiers may differ with respect to other parameters, will be considered. 
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The considerations are limited at first to a case in which the two amplifiers have the same shape 
of gain function but with differing gain amplitudes. If H\{o)) is the complex voltage gain of the 
first amplifier, then the complex voltage gain of the second will be taken to have the form 



// 2 (o>, t)= [l + a(*)]#i(ft>). 



(38) 



where a{t) is a stochastic function. 

If x\{t) is the signal resulting from amplifying the signal from the source X (this signal is l(t)) 
with the first amplifier and x 2 (t) is the result of amplifying l(t) with the second amplifier, then a 
straightforward Fourier analysis shows that 



x 2 (t)=[l + a(t)]x 1 (t), 



(39) 



and a similar relationship occurs for the signals y x {t) and y 2 (t) that arise from the signal m(t) from 
the source Y. This notation is summarized in figure 2. 

The most general form of the autocorrelation function of the sum of two random functions is 



ifj a ±b=[a(t)±b(t)][a(t + T)±b(t + T)l 



= i//„ + \\} b ± a(t)b(t + r) + b(t)a(t -f r). 



(40) 



From this it follows that 



</>U,+ */,) U 2 -y 2 ) = ^1*2 + ^iy 2 + ^i/,j-2 + *l*uiU2 -^\Xijxjt 



From (39) it follows that 



-x x {t)y x {t^7)ylt)x 2 {t^T)^y x {t)xAt^T)x 2 {t)y2{t J rT). (41) 



tyx x x % = xi{t)x 2 (t)xi(t + T)x 2 (t + r), 



--x%t)x%t + T)[l + a(t)][l + a(t + T)], 



(42) 



By means of a similar line of reasoning, every term in (41) can be put into a form similar to 
that in (42). Thus, 



^kxt+yMxt-yj — 0U, + .*/i)(.ri-i/,)(l + *fa + 2a), 



(43) 



I + m + n 



v 1 = x 1 + y 1 + z 1 Smoothing 



v 2 =x 2 -y 2 + z 2 

I -m + n 2 



Filter 



K3 



FIGURE 2. Radiometer circuit with terminology for sec- 
tion 3. 
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again 



and in the same way 



^i(x 2 -y 2 ) = (^ 1 x ] + $z iyi )( I + </'« + 2a), (44) 



^z^x 2 -y 2 ) = (^z\x i + ^j//,)( • + </>« + 2a), (45) 



ty Zl z 2 = tyz lZ {\ + i// a + 2a), (46) 

where zj is the signal that would result if the noise signal n 2 (t) (the equivalent noise generated in 
the input of the second amplifier) were amplified by the first amplifier. 
When these terms are all collected together, the result is 

i// rir2 = (1 + i|fa + 2a){{x\- y?) 2 + ^ + ^ 4- 1//„„ + ^ iyi + ^ + i//„. (/l + i^} • (47) 

Notice that in this case the second and third terms of (41) are cancelled by the sixth term. 

In order to simplify the analysis, we will consider two cases: (1) a a constant, independent of 
time, and (2) a(t) a random variable with a = Q. 

Case 1. When a is a constant independent of time, it is elementary to show that 



a= a 



(48) 



Thus 

0„ Bj = (l + a) 2 {(*?-y?) 2 + 2(»J£ i + ^ ) + <fc - +*,„ + * + * +*,,,}• (49) 

1 * 1 ' '' ' ' 11 "" 1 " l 11 

Therefore, both the d-c and the fluctuation terms are multiplied by (1 + a) 2 . Since this affects 
both the fluctuations and the deflection by the same amount, the sensitivity is not changed. 

Case II. Assume that a(t) is a random function of time, such that a = 0. At balance the term 
(x i — yf) — 0, so using a prime to denote quantities evaluated at the balance point, the result is 

^ I „ = (l+*.){2(*i I + i^ 1 ) + *. 1 » 1 + * Ml + *,,, i + ^ liiri + * llil , i }. (50) 

From this the spectral density is 

<?r,r, =d + <?«) *{2(£,, * <?*, + <?„, * £,, ) +Q. l *(Qx l + Qw l ) + Q ti *(Qx l + Q» ,) + <?,.*<?*,}• (51) 

Notice that (51) is similar to the result that is obtained in the case of identical amplifiers; the only 
difference is that the present result is the convolution of the original result with (l-\-Q a ). 
The spectral density of the multiplier output is 

<?« = « 2 <?vv 
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where Qr v is given by (51). The spectral density at the output of the low pass filter is 
lT(f) = G(f)Q u (f), 

1 

(52) 

Following the same line of reasoning that was used from (11) to (14), the root mean square fluctua- 
tion on the output can be seen to be 



w„ 



fx "I 1/2 

J Wifidf , 



+Q\ l *(Qx 1 + Qv 1 )+Q z{ *Qz 1 }df 



1/2 



(53) 



To see the effect of the gain fluctuations, consider that the amplifiers have a square bandpass 
of bandwidth B. Then, 

Q Xl *Q x =2Q*R*B, 

where R 1 (a)) = H 1 (o)) tf* (o>), 

is the power gain of the amplifier, and 



Q**Q,i*Q* x =\~ 2Q>R>BQ a (f-f')df'=2QjR\Ba>. 

J — DC 



(54) 



Thus the term (1 + Q a )*(Q Xl *Qxd can be written as 2Q 2 R 2 B(l -f a 2 ). Each term in (53) will produce 
the same effect, so the result is 

Wrms=(l+a 2 ) 1/2 Wrms(0), (55) 

where w rms (0) is the fluctuation level that occurs when there are no gain fluctuations, i.e., when 
a = 0. 

Since gain fluctuations increase the fluctuation level of the output but have no effect on the 
average deflection, it is evident that this effect results in a reduction in radiometer sensitivity by 
the factor (1 +a 2 ) 112 . 

While the above expression for the reduction in radiometer sensitivity has been derived only 
in the case of white noise and amplifiers with square bandpass, examination of the integral that 
leads to (54) shows that the result will be qualitatively similar even in a less idealized case. 

The next consideration will be to calculate the sensitivity of a radiometer with differing ampli- 
fier gain functions (transfer functions). In general these amplifiers could have differing center 
frequencies, and also the shape of the gain-frequency functions could differ. 

The average square output can be expressed as 



«*«)«=[* G(f)Q„df (56) 
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where G(f) is the power transfer function of the low pass filter. The input to the filter is given by 

u = av\Vi = a(xi + j\ + Zi)(x 2 — y 2 + 22). (2) 

By means of the theorem expressed in (4), the autocorrelation function of u can be put into the form 

a- 2 \\) u = (xxx 2 - yir 2 ) 2 + (i//.r, + ^ + ifc 1 )(<fcr 1 + */>*/ 2 + <fe 2 ) 



+ *,(0* 2 (* + t) ^ 2 (0^(^ + T) + yi(0r2^ + T) y 2 (0yiU + r) 



-* 2 W*i(* + t) yi(%2(* + r)-*i(*)* 2 (*+T) y 2 (0yi(^ + r). (57) 
Recalling that 

Q x =RiQi=H x H*Q u 

and similar relationships exist for the other variables, the spectral density of u becomes 

+ Re{[(//i#*<M*(// 2 //?^ 

-[(H^HiQiMHiHtQm^^HtHtQiHHtHtQm)^. (58) 

The form of the last term arises from the theory of cross-spectra [Goodman, 1957; Korn and Korn, 
1961]. 

The complete expression for the spectral density of the multiplier output is obtained by 
evaluating (58). From the resulting expression the fluctuation amplitude of the radiometer can 
be obtained by following the same procedure used previously. 

Next it is necessary to calculate the average deflection resulting from a particular combination 
of input signals. This can be obtained by evaluating (2) of the simple theory in this more general 
case. 

First, consider 



u = av\V2 = ci(xi-{- j\ + Z\)(x2 — y 2 + z 2 ). (2) 

Of the six random variables that appear in the right-hand term of (2), the only correlations 
that are not zero are the partial correlations of x\ with x 2 and of yi with y 2 . Therefore this expres- 
sion reduces to 

u = a(xix 2 — yiy 2 ). (59) 

Equation (59) must be evaluated in terms of the transfer functions of the two amplifiers, 
//i(/)4-// 2 (/), and the spectral densities of the noise outputs of the two noise sources. The Fourier 
transform of the random function / is 



S l {f)=[l{t)e- 2 «mt. (60) 

In the theory of cross-correlations, it is shown that [Korn and Korn, 1961] the time average of 
the cross correlation of x\ and x 2 can be expressed as 

<*,*>« f Re Lim i {Sxi S$J df= Lim ^j & x ^S t +SSJS^=\J(fI l Hf^HfH,)Q t 4f, (61) 
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where the symbol ( ) is used to indicate a time average. The final form of (61) results from the 
fact that S x = H\Si, etc. Also, due to the fact that the random processes are ergodic, the time 
average in (61) may be replaced by an ensemble average. 
In the same way 

^ 1 = ^(H,H* + H- z H*)Q m df. ■ (62) 

Therefore, 

u = a[^^ 2 -^ z ]^j(H l H^H 2 Hf)(Q l -Q m )df. (63) 

The above expression leads to the following conclusions: (1) In case the two sources have the 
same form of spectral density over the amplifier bandwidth, balance occurs independent of the 
amplifier response curves. Thus, for white noise, balance occurs when Ti = T m . (2) In case Qi 
and Q m have different forms, balance occurs where 

J (#,//* + H t H*)Q l df=j(H 1 HI + H- z H*)Q m df. (6 4) 

An expression for Au, the change in u due to a deviation from balance, is required. In the 
case where both signals consist of white noise, 

Q,-Q m = \k{T,-T m ). (65 ) 

Next, the assumption is made that 

T m =Ti-AT, 

then, Q,-Q m = \k[T,-(T,-AT)\ = \kAT. (66) 

In this case 

Au = |M7 , f(// 1 //* + W 2 //*) ( //. (67) 

If the smoothing filter has a response function at zero frequency given by Go, then 

Aw = GfAH = - GfkAT I (//,//* + Htf *)df. (68) 

8 j 

Again the assumption is made that the minimum detectable change in deflection is a change 
equal to the rms value of the fluctuations. The minimum detectable value of A7 1 is evaluated by 
equating the right-hand side of (68) and the square root of w 2 . Using the fact that 



1/2 

Wrms = (w 2 ) 



/: 



W(f)df 



the result is 



| X G(f)Q' u (f)df\ m = ^ GJ/W j (Htf* + H 2 Ht)df. (69) 
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When this solved for AF, the result is 

f* G(j')Q> u (f)df\ 
AT- 



^G^kl (H x H* + H 2 H*)df 



(70) 



Equation (70) is the general expression for the sensitivity of a correlation radiometer with 
amplifiers with differing complex gain functions. In order to interpret this result and obtain a 
feeling for the importance of the various parameters, it is convenient to look at two examples. 
First the case of two amplifiers with square gain functions centered on the same frequency fo; the 
first amplifier with bandwidth B\ and gain R \ and the second with bandwidth B 2 and gain R? 2 will be 
considered. To be specific, it will be assumed that B\ > B>. It will be assumed that no phase 
shifts occur in these amplifiers. 

In order to evaluate (70) it is necessary to evaluate the various terms of (58), recalling that the 
numerator of (70) contains only fluctuation terms; i.e., any d-c contributions must be ignored. It 
is reasonable to evaluate the convolution integrals that occur with the restriction that /is small 
enough that B> is enclosed within B\. When these evaluations are carried out with the assumptions 
that the radiometer is balanced and also that the low pass filter has a square passband of width b, 
the result is 

[numerator of (70)]= aG^b^BiR^R^kiTi^Tn). 

When this is compared to the similar result from the simple theory, for example (14), it is evident 
that the bandwidth that sets the fluctuation level is B 2 , the narrower of the two amplifier bandwidths. 
Also, the amplifier gains appear as the product R^R". 

The denominator of (70) must also be evaluated under the same conditions. This involves 
evaluating the integral 



/ 



{HxHf + HiHDdf, (71) 



where H\ and H 2 are the discontinuous real functions as defined above. When this is carried 
out, the result is 

[denominator of (70)] = £ GJ /2 k(4R° 1 R° 2 B 2 ). 

o 
These results lead to 

AT=2(£j U (Tt+T,,). (72) 

Thus, in this case the same sort of expression is obtained as in the case of the simple theory 
except that the bandwidth of the narrowest amplifier sets the sensitivity. 

Finally the case of two amplifiers with square bandpass with identical bandwidth but tuned to 
slightly different center frequencies will be considered. Again phase shifts will be neglected, the 
complex gains will be treated as real positive functions, and the radiometer will be assumed to be 
balanced. The resulting sensitivity is 

&T=2UjX\T l + T n ) (73) 

where B {) is the "overlap bandwidth" or the frequency interval common to both amplifiers. 
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4. Effect of Variation in Time Delay and Phase in the 
Amplifiers of a Correlation Type Radiometer 

Here it will be assumed that a differential variation in time delay can occur to the signals in 
the two amplifiers. Thus, at any instant the multiplier is comparing voltages corresponding to 
two different instants of time. Therefore, the multiplier output can be expressed as 

u(t) = v t (t)v 2 (t + T), (74) 

where r is the time delay difference in the two signal channels. In order to simplify this analysis, 
the assumption will be made that the two radiometer channels are identical except for the dif- 
ference in time delay. Thus the voltages that are multiplied together are 

vi(t) = x(t) + y(t) + zi(t), 

v 2 (t + T) = x(t + T)-y(t + T) + Z2(t + T). (75) 

Equation (7) shows that the fluctuation level of the radiometer depends only on convolutions 
of power spectra. Since these are not a function of r, the fluctuation level does not vary as a result 
of a differential time delay. Therefore, the change in sensitivity will be related to the change in 
average deflection caused by the time delay. 

The change in average deflection can be obtained from 



u(t, r) = [.xU) + yU) + 2iU)] [x(t + T)-y(t-{-T)-i-Z2(t+T)] 



= x(t)x(t + t) - y(t)y(t+T) = xjjAt) - Wt). (76) 

The Wiener-Khintchine theorem allows this to be put into the form 



u{t, r) 



j(Q,-Qy)e^df (77) 



In order to be able to evaluate the above integral, it is convenient to consider that both sources 
generate white noise and the amplifiers have square bandpass of width B centered on f {) . Thus, 
in this case, 



Ti(T) = (Q. r -Q u ) 



~ f-/ +fl/2 [L+BI2 

e iQ)T df+ e i(tiT df 

J-fo-B/2 Jfo-B/2 



^ n, ^ /~\ x t -i r sin tttB .„' 

= 2B(Q X - Q y ) cos 2ttt/o fl ■ (78) 

This can be put in more convenient form by noting that 

~u(0) = 2B(Q J -Q y ). 

Thus 

- - sin tttB 

u(t) = u(\j) cos Ztttjq — 

tttB (79) 

For the usual case in which B <^ f () and r is on the order of l//o, this can be approximated by 

u(t) « u(0) cos 2ttt/o. (80) 
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Thus, in the case of white noise sources and square bandpass amplifiers, time delay variations 
will not influence the balance point. They will affect, sensitivity, however. For maximum sen- 
sitivity, it is important to keep the difference in time delay in the two amplifiers small enough 
that cos 2ttt/o is approximately unity. 

A somewhat related question is the effect of a differential phase change in the two amplifiers. 
Again assuming white noise sources and amplifiers that arc Identical excepl for the phase shift, 
the average deflection can be computed by means of (63). With the assumption thai 



R 2 = Rie 



/</> 



(81) 



this equation becomes 



ui&^-iQt-QrrdRlV COS (f>). 



This can again be expressed in the form 



u((f)) = u(0) cos (f). 



(82) 



5. Some Effects of Using an Imperfect Multiplier 

The simplest form of multiplier to multiply two microwave signals together is probably that 
shown in figure 3. The input signals (v\ and v-i) are applied to two opposite arms of a hybrid junc- 
tion. The other pair of opposite arms are terminated by square law detectors. The outputs from 
these detectors are applied to a difference amplifier. In case the hybrid Tis matched, the detectors 
have the same sensitivity and their reflection coefficient is zero; and, if they are assumed to have 
a perfect square law envelope response, the output voltage z will be 



( V1 + V2Y ( v\ -v>i\ 2 



(83) 



where a is the coefficient of proportionality in the multiplier law. Thus, 



2 = 2av\ v-2 . 



(84) 



In an actual multiplier constructed on this principle, the ideal conditions suggested above will 
be only approximated. Thus, the output will not be an accurate multiplication of the input signals. 
A scattering matrix analysis of the multiplier junction indicates that if the junction is matched and 
perfect, the wave arriving at the detectors is equal to that assumed in the ideal case as long as 
either pair of leads are terminated in matched loads. Since the detectors (if bolometers) can be 
quite accurately matched, this will probably not be the phenomenon that sets the multiplier 
accuracy. This accuracy can depend on the accuracy of the "square law" of the detectors, 
however. 



v^ 



Jl 



Figure 3. Microwave multiplier. 



\ 



Differential 
Amplifier 
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In order to investigate the effect of the detector law on the multiplier accuracy, assume that 
both detectors follow the same law. If one detector is more sensitive than the other, the incoming 
signal can be attenuated, so that, referred to the input of the attenuator, the two detectors are 
assumed to have the same law with the same numerical coefficients. Thus, the voltages from the 
detectors can be expressed as power series 






(85a) 



h= ^ an Vvf. 



(85b) 



Then the multiplier output is 



, , ^ r ( vi+v£ \" ( V1-V2 Y 

These terms may be expanded with the binomial theorem; the result is 



(86) 



21/2 



2 2/2 



V2 + ^ (2vii*) + ^ (3v 2 v 2 - 



g3 
2 3/2 



a 4 



4) + gf(4»f% + 4i;itD + 



(87) 



The smoothing filter provides an estimate of the average of the multiplier output; thus it is neces- 
sary to compute z. This is, 



z = 2[a 2 {viV 2 ) + #4(^2 + V\V§ -f . . . ] , 



2[a 2 (x + y + z x )(x- y+ z 2 ) + a 4 {{x + y+ Zl )\x- y + z 2 ) + (x + y + Zl ){x- y+ z 2 ) A } + . . .], 



= 2[a 2 (x 2 -y 2 ) + 2a 4 (x 4 -y 4 ) + 3a 4 (x 2 -f)(z 2 l + z 2 2 )+. . .]. 



(88) 



The result is that ii x 4 — y 4 — at the same conditions that x 2 — y 2 — 0, the balance condition 
is not affected by any term up to the fourth in the crystal law expansion. This would be the case 
if both x and y are signals with the same statistical properties; for example, if both are thermal 
noise sources, they both possess Gaussian statistics. 

However, in the case that the two signals have different statistical properties, (x 4 — y 4 ) would 
not necessarily equal zero when (3c 2 — y 2 ) does. Thus, in this case, it is important to select de- 
tectors that are accurately square law. 

6. Conclusions 

The sensitivity of a correlation radiometer has been computed under a variety of conditions. 
In section 2 a very simplified calculation is used to derive the usual expression for the sensitivity 
of the radiometer. This is followed by a calculation of the sensitivity of a radiometer of the type 
suggested by Alfred (1962) in which the unknown noise signal is balanced against a sinusoidal 
reference signal. By performing these two calculations in an analogous manner, it is particularly 
easy to compare the sensitivities that result. In the case that the output fluctuations are pre- 
dominately due to the input noise signals, Alfred's radiometer is more sensitive than the conven- 
tional circuit that compares two noise signals. 
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The first part of section 3 demonstrates that it is not necessary for the two amplifiers to have 
the same gain. The calculation shows that the product of the gains is the parameter that deter- 
mines the output amplitude, instead of the individual amplifier gains. Next, it is shown that gain 
fluctuations produce a decrease in sensitivity. The last part of this section is concerned with the 
effect of using two dissimilar amplifiers. A general expression is obtained for the sensitivity of a 
radiometer with amplifiers with arbitrary gain functions. This expression is evaluated for the case 
of amplifiers with square bandpass and no phase shifts. It is shown that under these conditions, 
whether the two amplifiers have the same bandpass and different center frequencies or whether 
they have the same center frequency and different bandwidths, the sensitivity is determined by 
the "overlap" bandwidth. 

Section 4 deals with the effect of a differential time delay or differential phase shift in the two 
radiometer channels. It is shown that both of these effects result in a decrease of sensitivity. 

Finally in section 5 it is shown that, if the multiplier carries out the multiplication operation 
by forming the difference of the squares of the sum and difference of the two input signals, errors 
can result if the "square law" elements do not have a perfect square law response and if the two 
input voltages have differing statistical distributions. Thus, the radiometer proposed by Alfred, 
in which the comparison signal is sinusoidal, places a more stringent requirement on the multi- 
plier than does the more usual correlation radiometer, in which both input signals are Gaussian 
noise. 
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8. Appendix- Evaluation of Convolution Integrals 

Assume a random signal with white spectral density amplified by an amplifier with square 
passband as shown in figure A-l. 

If the spectral density of the input signal is Q as measured, the amplitude of Q x is RQo/2, 
where R is the power gain of the amplifier. The division by 2 occurs because the spectral density 
is assumed to be split equally between the positive and negative frequency regions. Then, 



Qx*Q; 



= \ Q * 



(f-f)Qx{f')df. 







Figure A-l. Square spectral density. 



FIGURE A-2. Functions appearing in integrand of con- 
volution integral'. 
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The value of this integral for small values of the parameter /is required. The term Qx{f~f) has 
the same shape as Q x (f') except that it is displaced an amount / and inverted on the frequency 
axis, as shown in figure A— 2. Since the integrand is the product of these two functions, the value 
of the integral is the product of the amplitudes of these functions times the frequency interval 
over which the product is nonzero. Thus, 

Q x *Q x = (R^j2{B-f) (A-l) 

for f<B. 

In this paper, the value of this integral is required under the condition that f*£ B. To this 
approximation, the result is 

Q x *Q x = 2b(rQ£J- (A-2) 

In general, the convolutions of the other spectral densities appearing in this paper are evalu- 
ated in the same way. The reason that the various convolutions of power spectra that are used 
have differing numerical coefficients is that the hybrid junction at the input of the radiometer 
divides the input powers; however, the noise powers appearing due to the effective input tempera- 
tures of the amplifiers do not undergo this power division. 

(Paper 71 C 1-253) 
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